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FOREWORD 


This  report  stems  from  a study  performed  spe- 
cifically for  Gibbs  and  Cox,  Inc.,  under  an  Office 
of  Naval  Research  sub-contract,  in  wnich  the  stability 
of  a particular  configuration  previously  tested  at  the 
Experimental  Towing  Tank  was  analyzed  (Reference  1). 
The  present  development,  which  was  performed  under  ONR 
Contract  No.  Nonr  263,  Task  Order  01,  is  more  general 
and  detailed,  and  takes  into  consideration  an  addi- 
tional degree  of  freedom  (surging  motion). 
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SUMMAftf 


The  equations  of  motion  describing  the  longitudinal  stability  char- 
acteristics of  an  uncontrolled  tandem  hydrofoil  system  in  smooth  water  with 
three  degrees  of  freedom  are  developed  in  the  present  report.  The  stability 
derivatives  in  the  equations  are  evaluated  on  the  basis  of  theoretical  hy- 
drodynamic formulas,  thereby  eliminating  the  necessity  of  their  experimental 
determination.  Comparison  of  the  motions  obtained  from  solution  of  the  equa- 
tions with  those  determined  from  previous  experimental  stability  tests  in- 
dicates a duplication  of  the  motion  experienced  in  the  tank  tests.  Varia- 
tions in  the  location  of  the  longitudinal  center-of-gravity  position  are 
shown  to  lead  to  different  types  of  motion,  which  is  in  agreement  with  the 
results  of  standard  aircraft  stability  analyses. 

The  additional  degree  of  freedom  in  longitudinal  surging  is  found  to 
indicate  a decrease  in  stability.  This  mode  of  motion  is  of  great  importance, 
however,  for  a study  of  stability  in  waves  where  the  waves  cause  variations 
in  the  forward  speed. 

The  equations  developed  in  this  report  may  be  extended  to  studies  of 
motion  in  waves  and  controlled  motion  by  incorporating  the  forces  and  mo- 
ments due  to  wave  motion  and  the  dynamics  of  the  control  system  in  the  equa- 
tions. By  assuming  the  resultant  system  of  equations  to  be  linear,  the  analy- 
sis may  be  performed  by  using  the  Laplace  transform  technique,  in  the  same 
manner  as  was  done  in  the  present  study. 
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INTRODUCTION 

The  use  of  hydrofoils  fcr  the  support  of  craft  operating  on  the  sur- 
face of  the  water  appears  to  be  the  most  promising  way  of  attaining  higher 
speeds  without  prodigious  amounts  of  power  or  without  a prohibitive  be- 
havior in  a seaway.  It  is  for  this  reason  that  the  hydrofoil  principle, 
while  not  new,  has  received  considerable  attention  recently. 

A hydrofoil  craft  in  the  foil-borne  condition  is  similar  to  an  air- 
plane in  flight.  Therefore,  it  is  not  surprising  that  the  longitudinal 
stability  of  hydrofoil  boats  is  one  of  the  problems  confronting  the  de- 
signer. A certain  amount  of  work  has  already  been  published  on  the  sta- 
bility of  motions  of  hydrofoil  systems  in  smooth  water  (see  References  2, 

3,  and  ii),  with  the  hydrodynamic  derivatives  evaluated  from  experimental 
data  and  the  downwash  of  the  forward  foil  of  a tandem  system  obtained  flom 
the  aerodynamic  theory  of  a foil  far  removed  from  the  surface.  Subsequent 
studies  of  the  hydrodynamics  of  foils  operating  near  the  free  surface  have 
yielded  theoretical  expressions  for  the  lift  and  drag  forces  and  the  down- 
wash  which  are  close  approximations  to  the  actual  effects  determined  ex- 
perimentally (References  5,  6,  7).  TTie  downwash  at  various  distances  aft 
of  the  forward  hydrofoil  may  even  be  an  upwash,  depending  on  the  forward 
speed  and  lift  coefficient.  Alterations  in  the  downwash  can.  of  course, 
a&terially  affect  the  longitudinal  stability  of  hydrofoil  craft. 

The  present  report  deals  theoretically  with  the  longitudinal  sta- 
bility of  fully  submerged  tandem  hydrofoil  systems  in  smooth  water  with 
three  degrees  of  freedom,  and  utilises  thecretical  expressions  for  the  hy- 
drodynamic derivatives.  It  is  a first  approach  to  the  problem  of  developing 
methods  of  theoretical  treatment  for  the  more  general  case  of  longitudinal 
stability  of  hydrofoil  systems  in  rough  water.  The  notation  used  has  been 
adapted  from  that  presented  in  Reference  d and  affords  a simple  treatment 
for  this  complex  stability  problem. 

In  addition,  the  results  of  some  previous  experimental  investigations 
are  compared  with  those  obtained  from  the  present  theoretical  study. 
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A = aspect  ratio 
C “ propeller  "drag"  coefficient 

t\0 

Cj.  * D/~  S'v  ^ , drag  force  coefficient  in  equilibrium  condition 

C.  = L/^  SV  2 , lift  force  coefficient  in  equilibrium  condition 

c.  G 

Cm  = SV  ^c  , pitching  moment  coefficient 

Cf  - T/^  SVg^  , thrust  coefficient 

c = mean  hydrofoil  chord 
D =•  drag 

d * vertical  distance  of  foil  to  reference  axis  of  hydrocraft 
g - acceleration  due  to  gravity 

H = vertical  displacement  of  hydrofoil  from  equilibrium,  positive  downward 

h = submergence  of  foil  below  smooth  water  surface 
2 2 

ig  " Kg  /X  » moment  of  inertia  coefficient 

Kg  = radius  of  gyration  of  hydrocraft  about  the  lateral  axis 

k m coefficient  of  accession  to  inertia 
L « lift 

« distance  from  l/U-chord  point  of  front  foil  to  1/U-chord  point  of 
rear  foil,  measured  along  the  reference  axis 

M = pitching  moment 

m , m^,  etc.  * stability  derivatives 

q 3 angular  velocity  in  pitch  (radians/sec, ) 

q 3 qt  , nondimensional  angular  velocity 

S - total  hydrofoil  area 

S„  “ area  of  forward  foil 
r 

SD  3 area  of  real'  foil 
K 
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T = propulsive  thrust 


t = time,  sec. 


£ = W/gpSV  , unit  of  aerodynamic  time 


u = increment  of  velocity  along  the  x-axis  in  disturbed  flight 


u = u/V  , nondimensional  increment  of  velocity  along  the  x-axis  in  dis- 
turbed flight 


V = resultant  velocity  of  hydrocraft  in  disturbed  flight 


= velocity  of  hydrocraft  in  equilibrium  condition 


W = all-up  weight  of  hydro craft 


increment  of  velocity  along  the  z-axis  in  disturbed  flight 


w/Ve  , nondimensional  increment  of  velocity  along  the  z-axis  in  dis- 
turbed flight 


axis  fixed  in  hydrocraft  in  disturbed  flight  in  direction  of  motion 
in  equilibrium  condition 


xj^=  distance  along  reference  axis  from  C.G.  to  l/U-chord  point  of  front 
foil 


xu,  X , etc.  = stability  derivatives 


lateral  axis  fixed  in  hydrocraft 


z = axis  normal  to  x-y  plane,  directed  downward 


z , z . z , etc . 
u v’  q 


stability  derivatives 


hydrofoil  incidence  measured  from  zero  lift 


-H/Vet  , nondimensional  height  increment 


small  increment 


downwash  angle  at  rear  foil 


angle  of  rotation  of  the  hydrocraft  in  pitch  from  equilibrium  con- 
dition 


root  of  the  stability  quartic  (or  quintic) 
W/gpSJL  , nondimensional  mass  factor 


water  density 

t/t  , nondimensional  time 
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Subscripts 

e denoting  equilibrium  value 
F denoting  front  foil 
ft  denoting  rear  foil 
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BASIC  ASSUMPTIONS  UNDERLYING  THo  THEORY 


This  particular  study  is  limited  to  the  motions  of  a noncavitating 
tandem  hydrofoil  configuration  in  the  plane  of  symmetry.  The  hydrofoils 
have  equal  chords,  but  the  spans  may  differ,  allowing  different  areas  and 
consequently  different  load  distributions  on  each  foil.  The  water  surface 
is  assumed  to  be  smooth  except  for  disturbances  due  to  the  motion  of  the 
foils  that  may  alter  the  surface  conditions  in  the  immediate  neighborhood 
of  the  foils.  Time  lag  terms  are  utilized  to  express  the  influence  of  flow 
phenomena  occurring  at  the  forward  foil  upon  the  rear  foil  at  a later  time 
(due  to  moving  through  a distance  Jo). 

There  are  five  variables:  u , w , 6 , q , and  B . Thus,  in  addition 

to  the  three  usual  equations  of  aerodynamic  longitudinal  stability,  there 
must  be  two  other  equations  which  will  be  furnished  by  purely  kinematic 
cons ider  at ions . 

The  hydrodynamic  force  and  moment  coefficients  depend  upon  only  the 
nondimens ionalized  vertical  and  angular  displacements  and  linear  and  angu- 
lar velocities.  Terms  due  to  buoyancy  are  neglected  since  they  are  not  sig- 
nificant in  comparison  with  the  dynamic  forces. 

Tiie  hydrodynamic  lift  and  drag  force  coefficients,  CT  and  (i,  are 

D 

considered  as  functions  of  the  incidence  angle,  a,  the  speed,  V,  and  also 
of  the  depth  below  the  water  surface.  These  coefficients  are  to  be  evalu- 
ated at  the  equilibrium  position. 

The  effect  of  thrust  variation  with  speed  and  depth  can  be  taken  into 
account,  but  is  considered  to  be  negligible  in  the  present  development.  The 
elastic  distortion  and  slipstream  effects  on  the  rear  foil  are  to  be  ne- 
glected. The  thrust  is  assumed  to  act  through  the  center  of  gravity  and 
along  the  x-axis. 

In  addition,  all  the  simplifications  applicable  to  the  equations  by 
use  of  linearization  by  the  method  of  small  disturbances  are  inherent  in 
the  present  development. 
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AXIS  AND  COORD INAi.iS 


The  origir.  is  at  the  center  of  gravity  of  the  hvdrucraft,  The  x-axis 
is  forward  along  tne  water  direction  in  the  equilibrium  condition  and  is 
fixed  in  the  hydrocraft  during  the  disturbed  motion.  The  z-axis  is  down- 
ward in  the  plane  of  symmetry  of  the  hydrocraft  and  perpendicular  to  the 
x-axis,  while^the  y-axis  is  to  starboard.  The  positive  direction  of  the 
angle  of  rotation  in  pitch,  9,  is  counter-clockwise.  The  aocve  sign  con- 
vention, as  well  as  a representation  of  Lbe  forces  and  Velocities  asso- 
ciated with  the  disturbed  motion  of  the  hydrocraft,  is  shown  in  the  fol- 
lowing sketch:  __ 


HI"  I7.CWT> 


v ' OK  i i.  i_a’«  i i.  /i 
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EQUILIBRIUM  CONDITIONS 

When  the  hydrofoil  system  is  in  static  equilibrium,  it  is  moving 
forward  at  a velocity  Vg  in  the  x-direction  at  a definite  depth  of  sub- 
mergence of  the  foils  and  with  the  reference  axis  at  zero  trim.  To  achieve 
these  conditions,  the  following  requirements  must  be  satisfied: 

(a)  the  sum  of  the  forces  in  the  x-  and  z-directions  must  equal 
zero. 

(b)  the  moment  about  the  center  of  gravity  due  to  the  forces  on 
the  foils  must  equal  zero. 


These  conditions  are  expressed  analytically  as 
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EQUATIONS  OF  MOTION  IN  DISTURBED  FLIGHT 


Force  Components  Along  the  x-Axis 

The  forces  in  the  x-direction  in  disturbed  flight  are  related  by 
the  following  equation: 

| (5?  + v<^j  “ <T  + AT>  t (LF  * ^F  + LR  * 

- (0^  ♦ ADR  * DR  *ADR)cosw  -W  sin  0 . (3) 

Taking  into  account  the  equilibrium  conditions  (equation  (1))  and  ne- 
glecting small  quantities  of  higher  order  gives 

| ^ - AT  ♦ (Lp  ♦ Lr)w  - (ADp  ♦ ADr)  - W0  . (li) 

2 

Dividing  by  pSV  in  order  to  pass  to  nondimensional  form,  and  simpli- 
6 

fying,  yields 

dS  _ at  . 1 „ A adf  *adr  1 c 0 f?) 

— c w n 2 L8  ' 

pSV  ^ L pSV  1 L 

© © 

From  consideration  of  the  aerodynamics  of  propeller  thrust  (Reference  ti), 
the  following  equation  is  obtained: 

AT  . _c  A ^ » <*T  (6) 

TTT  AS  AS  pSV  dv  ' 


However,  it  will  be  assumed  for  the  present  study  that 


0 


The  increment  of  drag,  Ad,  is  due  to  the  perturbation  quantities 
and  is  found  from  the  following  considerations: 


S3  " | S(V,  ‘ u ♦ q*)2^  * 


ac  acn  ac  ac  \ _ 

^w  + aTe*3FH+3Tq/2 


2 

sv  4c: 

e D 
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AD  =l(c  +^DA+^DG+^H  + afDq)(1  + 2A  + 22l)  -Jr  . 

2\  D da  5 6 dH  dq  7 \ V / 2 D 


If  it  is  assumed  that  the  pitching  velocity  motion  causes  changes 
only  in  the  local  angle  of  attack  at  each  foil  and  does  not  influence  the 
resultant  velocity  at  the  foil,  then 


ADF  tADR  _ SF 


/V  '{y)f  '[-d)e°  ’Ur)," 


4h 

•*ffl  .•*)}•#©  .4®.V 


q • (?) 


The  dependence  of  drag  forces  on  q and  w are  of  insignificant  magni- 
tude and  hence  terms  expressing  any  relations  of  this  type  will  be  ne- 
glected. The  effect  of  downwash  on  the  rear  foil  and  its  influence  on  the 
components  making  up  the  stability  derivatives  is  discussed  below. 


The  contribution  to  the  drag  force  due  to  downward  velocity  w has 
the  same  effect  as  that  due  to  changing  the  angle  of  incidence  at  the  foals 
by  the  additional  angle  w = w/Vg  . Now,  taking  into  account  the  change 
due  to  w when  the  angle  of  incidence  at  the  rear  foil  is  expressed  as 


t w d<  / ^ w \ 

^"““h  Ve’da\eF  Ve) 


with  € = downwash  angle,  and  considering  changes  due  to  w only  gives 


w d*  / w 


R V da 
e 


i-£’ 

da 
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Thus  tne  component  expressing  the  dependence  of  the  increment  in  the  drag 
force  on  w is  given  now  as 


The  derivative  dC^/dB  will  not  be  determined  directly  as  such, 
since,  as  mentioned  previously,  the  hydrodynamic  derivatives  are  assumed 
to  be  functions  of  angle  of  attack,  speed,  and  depth  of  submergence.  This 
derivative  will  be  transformed  so  that  it  will  be  found  in  terms  of  de- 
pendence on  depth.  The  relations  among  the  derivatives  which  allow  G 
terms  to  be  written  in  terms  of  change  in  depth,  H,  and  subsequently  the 
r.ondimensional  height  increment,  8,  are  given  below: 


*CD 

w 


aCD  dH 

<3h  Je 


^DdB 

dp  de 


8- 


Frcm  Sketch  1,  -H  =*  x6,  and  by  matching  the  value  of  H frcra  the  defi- 
nition of  8 , it  is  seen  that 


V t8  - x8 
e 


e 


In  order  to  maintain  the  proper  orientation  of  the  axes  and  coordinates, 
x « x for  the  forward  foil  and  x « -(£-  xj^)  for  the  rear  foil  (8 
measured  positive  for  a counter-clockwise  rotation).  The  component  ex- 
pressing dependence  of  drag  force  on  0 is  then  given  ac 


1 


rs4*V) 

~\Tb) 


sr(*cd)  + 
T \W)  - “ 


- Xx 


v t 
e 


(10) 


The  depenc.ence  of  the  drag  force  on  the  vertical  displacement  H 
is  converted  to  a nondimens ional  dependence  on  8 . Thus  the  component 
of  drag  varying  with  height  is  given  as 
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Hydrodynamic  forces  due  to  a pitching  velocity,  q,  arise  from  the 
change  of  angle  of  attack  due  to  the  resulting  curvature  of  the  stream- 
lines. Derivatives  with  respect  to  q can  be  transformed  to  derivatives’ 
with  respect  to  qx/Vg  » or  in  effect,  with  respect  to  an  angle-of-attack 
change  due  to  the  pitching  velocity  of  a foil  at  a distance  x from  the 
center  of  gravity: 


n W 


dc„ 

_ D 

d&  VD 


Since  the  front  foil  has  an  upward  velocity  equal  to  qx^  , it  has  an  ap- 
parent decrease  in  angle  of  attack  of  qx^/Vg  . The  rear  foil  experiences 
an  apparent  increase  in  angle  of  attack  equal  to  q( £ - x.  )/V  » Thus  the 
changes  of  angle  of  attack  due  to  pitching  velocity  are  -qx^/V  for  the 
forward  foil  and  q(^.-  x^)/Vg  for  the  rear  foil. 


Because  of  the  time  lag  for  the  downwash  from  the  front  foil  to 
reach  the  rear  foil,  there  is  a downwash  angle  lag  at  the  rear  foil  due  to 
pitching  velocity  (Reference  3).  The  flow  at  the  rear  foil  at  time  t is 
influenced  by  the  downwash  which  is  created  at  the  front  foil  at  time 
t ~-^/Ve  . This  effect  is  in  addition  to  the  local  change  in  angle  of  at- 


tack at  the  front  due  to  pitching, -qx^/V 


Tnus , 


dc  [ 
da 


\ Vc  J V * 
v e / e J 


and  at 


. a*  r qa'  ^r 

€ « — a„  + n 

L F e J 

The  angle-of-attack  change  at  the  rear  foil  due  to  pitching  velocity  is 
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R 


o.(  Ji-  x^) 

V 

e 


d€  K-  Xl) 

da  V 
L e 


a 


R 


xt>  /,  d«\ 

V V1  da  / 

a \ > 


(13) 


By  taking  into  account  the  local  ai  gle  change  due  to  pitching,  down- 
wash  effects,  and  modifications  due  to  the  time  lag,  the  resultant  com- 
ponent of  the  drag  force  dependent  on  q is  given  as 


\i  - 

(a£ 

s \J7J  F jl  s' 

iL'  d0' 

JL 

[\  ' 

After  the  nondimensional  factor  fi  *»  W/gpSX  is  introduced,  then 

qi/V  = (l//x)q  . Thus  the  contribution  of  pitching  velocity  to  the  drag 
0 

forces  may  now  be  written  in  terms  of  the  nonaimensional  variables  as 


Associated  with  any  accelerated  or  decelerated  motion  of  a body 
through  a fluid  is  the  virtual  mass  of  the  entrained  fluid.  For  a body 
such  as  a hydrofoil  at  small  angles  of  attack,  the  virtual  mass  due  to 
acceleration  in  the  x-direction  is  negligible,  aid  consequently  no  forces 
due  to  virtual  mass  are  included  in  the  equation  for  the  x-direction. 


By  combining  the  various  components  of  the  forces,  including  suitable 
modifications  due  to  downwash  effects,  the  resultant  equation  for  forces  in 
the  x-direction  becomes 


du 

dr 


X U - X 0 

u w 


xee  - V 


0 


(16) 


where 
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1 

2 


> (17) 


Force  Components  Along  the  z-Axla 

The  equation  relating  the  force  components  along  the  z-axis  may  be 
written  as 

I [5?  ‘ (Ve  + u)q]  " -(LF  * ALF  * LR  *ALR)c°S" 

- (Dp  ♦ AD_  ♦ D + AD  )sinO  ♦¥  cos  6 . (18) 

i t n K 

Taking  Into  aocount  the  equilibrium  conditions  (equation  (1))  and  neglecting 
mull  quantities  of  higher  order  gives 


i (If  - v.->)  ' -toLF  * flV  - (df  * DE 


\A 

)w 


(19) 


Dividing  equation  (19)  by  pSVe  in  order  to  pass  to  nondimensicnal  form, 
and  simplifying,  yields 

.A 


dw  * ALF*ALR  3.  a 


3t 


pSV 


T~  ' ? cdw  * 


(*o) 


By  performing  the  same  type  of  analysis  used  in  deriving  equation 
(7),  the  following  terms  are  obtained: 
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The  components  of  the  lift  forces  given  in  equation  (21)  do  not  contain 
% 

any  quantitative  measure  of  the  effects  of  dcwnvash,  and  thus  have  to  be 
corrected  to  toics  account  of  this  effect*  Alscv  some  of  the  component  terms 
have  to  be  modified  so  that  they  can  be  expressible  in  terms  of  nondiman- 
si  onal  variables. 

Applying  the  downwash  correction  to  the  component  dependent  on  w 
gives  an  expression  for  the  lift  force  variation  analagous  to  equation  (9)* 


The  component  of  lift  which  is  dependent  upon  6 is  expressed  below 
in  terms  of  the  nondimensional  derivatives  with  respect  to  0 and  is  equal 
to 


1 

2 


fr  (*V\  xl 


fs  Oh)  x'  ‘ * 
s Vae !R-pf 
0 — 


e 


(23) 
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The  component  which  is  dependent  upon  the  vertical  displacement.  H 
is  converted  to  a nondimensional  dependence  on  8 and  is  expressed  as 


Variation  of  the  lift  with  pitching  velocity  q is  due  to  the  dif- 
ferent local  angle-of-attack  changes  at  the  front  and  rear  foils,  the 
downwash  effect,  and  the  time  lag.  By  introducing  the  effects  of  these 
factors,  the  contribution  of  the  lift  due  to  q is  given  as 


Using  the  nondimensional  factor  ft  and  converting  to  the  nondimensional 
variable  q yields 


There  is  a substantial  dependence  of  the  lift  force  on  the  accelera- 
tion variables  w and  q that  yields  terms  whose  magnitudes  are  suffi- 
cient to  be  included  in  the  z-force  analysis.  The  two  accelerations  w 
and  q produce  forces  that  are  due  to  the  virtual  mass  associated  with 
the  accelerations.  These  virtual  mass  terms  are  noncirculatory  potential 
flows  that  do  not  produce  an  induced  vortex  wake.  Thus,  there  will  be  no 
effect  on  the  rear  foil  due  to  these  forces  acting  at  the  forward  foil. 

In  addition  to  it3  influence  cn  the  virtual  mass  terms,  the  accel- 
eration w also  causes  forces  on  the  rear  foil  due  to  a time  lag  factor. 
The  lift  at  the  rear  foil  at  time  t is  influenced  by  the  vortices  which 
were  cast  off  by  the  front  foil  at  a time  t - X/V  . Thus,  with  the 
angle  of  attack  at  the  front  foil  expressed  only  in  terms  of  its  varia- 
tion with  w , 


f » 
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da 
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V 


a 


(27) 


lhen  at  time  t = C , the  effect  of  v on  the  angle  of  attack  at  the 
rear  is 


a.. 


^ ["« (■  0] 


d«  wX 
35  ? 
ve 


(20) 


The  lift  force  at  the  rear  due  to  this  time  lag  effect  is  given  as 

r , ( Ps„2  p 2^  <iCL  ^ d<  wX 

L \JT/  “R  2 SRVe  2 SRVe  h dH  “?  * 


(29) 


L is  converted  into  nondimensicnal  form  by  dividing  equation  (29)  by 

pSV  2: 
e 


psv 


l*(dcA  *± 

2S  \7a’  / da 


R V 

e 


(30) 


Now  w can  be  changed  into  the  ncndimensional  form  of  a«)/dr  by  use  of 


the  definitions  of  t and  , yielding  the  relation 
v 2 . 

• dw  1 e dw 

V*dt^Jd?  * 


(31) 


Then 


sT7 


pbV 


1 SR  f3ClA  d€  dO 
2/i  S \3a  / R da  3t 


(32) 


From  this  expression,  a stability  derivative  can  be  defined* 
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Forces  due  to  virtual  mass  will  now  be  considered.  Because  c£  the 
similarity  between  an  elliptic  disc  and  the  planform  of  an  actual  hydro- 
foil, the  theoretical  calculations  made  for  an  elliptic  disc  can  be  used 
to  determine  the  virtual  mass  of  the  hydrofoil.  Ibis  apparent  mass  is 
equal  to  the  mass  of  an  ellipsoid  of  revolution  with  the  span  as  axis  and 
the  chord  as  diameter,  and  of  the  same  density  as  the  fluid.  The  volume 
is  (U/3) na^b  , where  b is  the  half  axis  and  a the  largest  radius  of 
the  ellipsoid.  The  inertia  factors  for  the  correct  aspect  ratio  are  ob- 
tained from  Reference  9*  Thus,  the  apparent  volume  K is  equal  to 
k(li/3)*a^b  , where  k is  the  coefficient  of  accession  to  inertia. 

It  is  possible  to  obtain  a foil  having  the  same  area  and  same  ap- 
parent volume  as  the  elliptic  disc,  but  it  will  have  a somewhat  different 
aspect  ratio.  The  apparent  volume  of  the  foil  is  then  K * k(n/l*)Sp  Rc  . 


Thus,  with  cE_ 


S * nab  (where  b - span  of  hydrofoil),  c = (l6/3n)a. 
F»« 


Since  the  inertia  factor  k has  to  be  the  same  in  both  cases,  it  depends 
upon  the  aspect  ratio  of  the  elliptic  disc  which  is  equal  to  ijb/na  . The 
aspect  ratio  of  the  elliptic  disc  is  found  to  be  1.166  times  the  aspect 
ratio  of  the  hydrofoil.  Since  the  inertia  factor  k does  not  vary  much 
for  aspect  ratios  greater  than  6.  and  since  the  factor  expressing  the 
ratio  of  the  aspect  ratios  of  the  disc  and  foil  is  not  much  larger  than 
1.0,  the  values  for  k given  in  Reference  9 for  the  disc  aspect  ratio 
may  be  used  without  any  correction. 

The  apparent  mass  is  then  pk(n/li)S,,  Dc  and  the  forces  are  the  mass 

* t H 

times  the  acceleration.  The  force  due  to  the  vertical  acceleration  w is 
a lift  force  -Z'  , where 


'Z'f,R  " pk  JT  SF,RW 


The  total  lift  is  then 


-Z'  - pk  I =<SF  . SB)i  -pkfsci.pkjsci^g 


t £ 7T  3?  • 
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Dividing  equation  (35)  by  pSVp  to  make  it  nondimens ion a]  yields 
stability  deriva^ ive 


2 w “ - 1 1 * (36) 

The  force  due  to  the  virtual  mass  which  is  contributed  by  q is 
considered  as  a lift  force,  and  is  expressed  as 

l'f,r  = p T SF,nc^xF,R  * (37^ 

This  force  acts  at  the  center  of  area  of  the  surface,  and  since  longi- 
tudinal reference  dimensions  are  measured  from  the  quarter-chord  point,  a 
term  must  be  added  tc  give  the  proper  distances  to  where  the  force  acts. 
Therefore, 


, . kn  _ • 

L p y SPC<? 


i («*-?)  *"t  •**(*- 


" p if  Sc  ' T (xi.  * f)  * T (* * xl  * k)  * * (j 

The  angular  acceleration  q may  now  be  expressed  in  terms  of  r.cndimen- 
sional  variables  by  the  relation 


dq  X 2* 

e 


TllU3, 


L'  - p 


bp2  Ji' 


J - T (XA  ' c)  * -£  (L~  xt  * c)_.5? 


L'  1 knc  l SF  Xl  " 4 . SR  XA*  u^aQ 

rr  " “7  HT  S 7 T X /3r 


Psve  M 


Another  stability  derivative  is  now  defined  as 
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The  variois  components  of  the  z-forces  will  now  be  combined,  with 
the  virtual  mass  contributions,  time  lag  factors,  and  dcwnwash  effects 
included,  yielding  the  final  relation 
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Moments  About  the  y-Axis 

The  moment  equation  for  pitching  about  the  v-axis  is  derived  below 

in  nondimensionai  form.  The  resultant  moment  is  -:sced  as  the  sum  of 

equilibrium  and  hydrcdvnanic  contributions : 

\ 


M = ft  ♦ M u + M w + + M,3  + M q ♦ M.w 

e u w 0 3 q v 


Equating  this  expression  to  the  product  of  moment  of  inertia  and  angular 

2 

acceleration,  expressed  as  (W/g)Kn  (dq/dt)  , and  subtracting  the  equi- 

D 

librium  terms  which  are  zero,  gives  the  following  equation: 


- K 2 ^ - M u - M w - M.G  - Hi  - H q - M.i  * 0 
g b dt  u w 0 q'1  *w 

2 2 2 

This  equation  is  then  multiplied  by  u.  /pSjttf  L , where  i =*  Kc  /// 

K being  the  radius  or  gyration  of  the  hydrocraft  about  the  y-axis: 
B 


dq  p Mu 

-f  - u - 


gpstv/  dt  ps/»e  iB 


pSiVe  :,B 


psxve  iB 


Psx»e  iB  Ps^ve  iB  Psive  iB 


A A A 

Now  with  q * qt  and  t * tr  , then 


aq  M t « JL  IS 
^ tdr  f2  dT 


and  with  t and  p defined  as 


W 

^ gpSl 
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the  following  rel.-:,i.on  results: 


. iL  _H Q . i£ 

. r l ir  u j 


r iBpSlVe 


iB  Psive 


B OS.4V. 


£.  P 3 


psz  veiB 


Equation  (U7)  may  now  be  written  at 


um  . 

’ u a “ w A 


u m • a 

Tf  C ■ 0 ’ <“> 


where  m , m , n„  , m„  , m , and  m*  are  nondimensional  partial  de- 
u w 0 0 q w 

rivatives  (stability  derivatives). 

The  terms  m , m.  , m.  , and  m may  be  easily  found  in  terms  of 
w ’ 6 0 ’ q 

the  moment  coefficient  . The  moment  is  put  into  nondimensional  form 

by  dividing  by  pSV  2 Jt  , giving  M/pSV  ^ . However,  C = M/(p/2)SV  2c 
© g m G 

and  the  two  forms  are  related  by  M/pSV  “U  « (c/2jt)C  . Thus, 

e m 

. dC 
_ c m 


mTld 


dC 

c m 

* -oT-  -TS- 
U c 4,  (fv 

dC 

c m 

m8 " tz  aF 


w 

aTl 


In  later  sections  of  the  report,  the  term  will  be  found  and  shown  to 
be  zero  due  to  the  equilibrium  conditions,  and  the  term  m^  due  to  time 
lag  will  also  be  determined.  In  addition,  moments  due  to  the  lift  forces 
generated  by  virtual  masses  will  be  found  and  included  in  the  final  moment 
equation. 


Security  Information 
CONFIDEN  UAL 


CONFIDENTIAL 
Security  Information 


KINEMAiiC  COrJLiJ  T.T  ; 


The  vvi.->  kiremotir  ° ■ ’■  ~«ce3s  • .u  order  to  have  five  equa- 
tion.0 . . ..  jr.  The  first  condition  expresses, 

tne  equality  c*  time  rate  of  change  of  pitch  angle; 


(ii9) 


It  is  now  necessary  to  express  the  relation  betvec:;  the-  vertical 
velocity  component  dH/dt  and  its  x and  2 components.  It  is  seen  from 
Sketch  1,  page  7,  that 

HU 

- — * (V  + u)sin6  - w ccs  6 , (50) 

dt  e 

or,  if  only  small  quantities  of  the  first  order  are  retained,  and  non- 
dimensional  variables  are  introduced,  then 

IS.  + 0 - 6 = 0 . (51) 
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LONGITUDINAL  STABILITY  DERIVATIVES 

In  this  section,  the  longitudinal  stability  derivatives  are  dism- 
issed in  detail.  The  various  hydrodynamic  derivatives  appearing  in  the 
expressions  for  the  stability  derivatives  will  be  evaluated  in  a later 
section  of  the  report.  Experimental  values  of  these  hydrodynamic  deriva- 
tives are  obtainable  from  towing  basin  tests. 

my  , Change  in  Pitching  Moment  with  Perturbation  Velocity 

The  moments  will  be  broken  down  into  drag  and  lift  forces  multi- 
plied by  their  respective  moment  arms. 

The  drag  and  lift  forces  due  to  the  longitudinal  perturbation  ve- 
locity u are  obtained  from  equations  (7)  and  (21),  respectively: 


Drag  Forces : pSV^C^u 


Lift  Forces:  pSV 

e L 


(52) 


These  forces  may  be  written  in  a slightly  different  form  in  order  to 
separate  the  actions  on  the  forward  and  rear  foils: 


Drag  Forces: 


Lift  Forces 


p°Ve2pFvK)5 


: pSV.2  C ♦ ^ c 'W 

c \s  lf  s lk/ 


> (53) 


The  moments  about  the  center  of  gravity  are  then  given  as 


M - pSV 


& v*  - T °> 


a-**)-  ^ - s’Vnj1 


2 [Vi  Vf  ■ “A] c 


(SU) 
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From  the  equilibrium  conditions  (equation  (2)),  M = 0 ; hence,  “ 0 

and  m = 0 . 
u 

m , Change  in  Pitching  Moment  with  ' .rtical  Velocity 

The  drag  and  lift  force?  Hue  to  vertical  velocity  are  found  from 
equations  (9),  (17),  (22),  and  (L3): 

Drag  Forces:  - § SV.2-^  ^ ‘ T (^)  J 


4* 


Lift  Forces:  ^ SV£ 


(55) 


fi  2«)pRc  "fljL 

V1  ' da  / [s  CDr  S \Jal  RJj  Ve 


The  moments  about  the  center  of  gravity  are  then  given  tv 


(i[cD/(vr)  ]xA-T  (x-S)  [Cdb* (*f)  **>}  f 

J'sjt  • (56> 


and  subsequently, 


-o- 


ac 

m 

3T 


* |{t[%  ' (jf)  Jvt  i1  ' li)  [%  • {77)  Jdn}  • (57) 
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But 


M = | SV 
w 2 e 


M*  I SVe{©} 


and 


M 


m 


w 


w pSXVr 


by  definition.  Hence, 


"»  ' Tf{©} 


(50) 


and 


mw  = “2X17 


dc 

c m 


me  , Change  in  Pitching  Moment  with  Pitch  Attitude 

A rotation  of  the  reference  axis  in  the  plane  of  symmetry  causes 
changes  in  the  depths  of  submergence  of  the  hydrofoils,  thus  altering 
the  drag  and  lift  forces  on  each  foil.  These  variations  can  then  be  con- 
verted into  a variation  in  the  pitching  moment  about  the  center  of  gravity. 

The  drag  and  lift  forces  due  to  submersion  changes  caused  by  rota- 
tion can  be  obtained  from  equations  (10),  (17),  and  (23): 


Brag  Forces:  tj  SVg2 


Lift  Forces:  | SV„2 


CL+S 


i(d%\  u V'M 

3 \w  / F v $ - s \w ) 

e 


X-  I } 

R 7% 


\{dS\  ?±-s-*(dcA  X~x7 

5 / F V % 3 ' R V \ 

— e e J 


The  term  in  the  drag  force  component  acts  through  the  center  of 
gravity  and  thus  does  not  affect  the  moments.  The  moment.?  about  the 
center  of  gravity  are  then 
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m„  , Change  in  Pitching  Moment  with  Vertical  Displacement  of  the  Center  of 
8 1 


Gravity 


Changes  in  the  drag  and  lift  forces  with  vertical  displacement  pro- 
duce moment  changes  about  the  center  of  gravity. 


The  drag  and  lift  forces  due  to  vertical  displacement  can  be  ob- 
tained from  equations  (ll)  and  (2lt): 


Drag  Forces 


e5.Esv2fe(V\  .*J b(V)1 

“s-  2 be  [s  \ap  ) F S V39/jJ- 


Lift  Forces;  & SV  2 


s \ds  l F s \as  / R 


3 . 


The  moments  about  the  center  of  gravity  are  then 


M - | SV,2 


[s  \HW/r?  **  S \wl 


U-  Xx)r  P 


&sv  2 


2 e 


fsF  d fa  (a°D^  A . 

{tVJsV/f  S Wi/jf  p 


(6) 

a .iD)  4 

C y S w / F F S \d  8 / R aRj 


Sr  ■ rf©V  if®l  • 


l J J 
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pSl'J 


by  definition.  Thus, 


R-h7$ 
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ac 

c m 

m9  = TL  7T 

, Change  in  Pitching  Moment  with  Pitching  Velocity 

The  pitching  velocity  causes  changes  in  the  angle  of  attack  at  each 
hydrofoil  which  in  turn  bring  about  changes  in  the  drag  and  lift  forces. 
These  force  changes  then  contribute  to  the  pitching  moment  about  the  center 
of  gravity. 

From  equations  ( ill ) and  (25),  the  drag  and  lift  forces  due  to  pitching 
velocity  arc  obtained; 


D 2 

Drag  Forces  £ SV 
2 e 


Lift  Forces  | SVg2 


(4(5),  M(S), [(■•£) 

(4(a,i-wa,&-s)^]L; 


Moments  about  the  center  of  gravity  are  then 


M = - SV  2 
2 e 


.■(-?[©,  i]vi  fa,  fa*  (-  a)]}4  if  <« 
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ana 


M 


m = 

q pS/fv 


■ft: 


e 


by  definition.  Thus, 

"o  - A(®|-  Ttf®}  <™> 

A J m.  Jr 

and 


m»  , Change  in  Pitching  Moment  with  Vertical  Acceleration 

The  importance  of  the  term  rrj  arises  from  the  fact  that  the  rear 
hydrofoil  at  time  t is  influenced  by  the  vortices  which  trere  cast  aft 


by  the  first  hydrofoil  at  time  t - -t/V 


Thus,  the  expression  for  lift 


on  the  rear  foil  given  in  equation  (29)  is  obtained: 
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The  moment  produced  by  the  rear  foil  is  given  as 
o ~ u 2 / de  wX  , i . 


KC  " ■ I V (■jrl  da 

H 


d€ 


(X.-  xx) 

By  definition,  fx  :.rj  • m^ps  jf  ; thus, 


m»  = - 

w 


1 

d« 

X‘  XJL 

2fL 

S 

x da: 

l 

y (7i) 


(72) 


Pitching  Moments  Due  to  Virtual  Mass 

The  accelerations  w and  q produce  forces  that  are  due  to  the 
virtual  masses  associated  with  the  accelerations.  Since  these  forces 
are  of  the  nature  of  lift  forces,  the  moments  due  to  such  forces  may  be 
determined. 

The  lift  force  due  to  the  virtual  mass  associated  with  the  accelera- 
tion w is  given  by  equation  (35): 


Lift 


-Z- 


P 


knc  ,r 

T wf 


♦ 3r)w 


The  moment  derived  from  this  force  is  found  by  use  of  a moment  arm  to 
the  center  of  area  of  the  foil  surface  and  hence  is  corrected  by  referring 
distances  to  the  quarter-chord  point  as  in  the  development  of  equation 
(38).  Moments  about  the  center  of  gravity  are  then 


M' 


knc 

p T" 


(x/  ‘ £ ) ‘ SR  (J 


s?lxi 


(73) 
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,,,  _ knc 

M' « = o 

w a 


s j (xf  ‘ t)  " T {£  " x^  + 


H-m'-  = — r 
pSX 


Thus* 


, _ knc  / Xl  ~ £ _ ^R  ^ yl  + 

1 w m\  s / s jl  / 


The  lift  due  to  the  virtual  mass  associated  with  the  angular  ac- 
celeration a is  given  by  equation  (3b): 


Lift  » L*  = p S 


Moments  about  the  center  of  gravity  are  then 


M!5pts;I  (x*  ' f)  ~ ~i  [f~ * i)  ] q 


M' • = - p 

q 


It  can  be  shown  that  m'^  = M'*/pS l?  , which  then  goes  into  the  moment 
equation  (Uo)  as  -(m* ^//iig  ) (dq/dr  ) . Thus, 


q 17  S 


kite  (Xf  U, 


SR  ‘ Xl  + 


The  final  moment  equation  including  the  effects  of  time  lag  and 
virtual  mass  contributions  is  then 

A - I'l)  4 a . f h e - 05. 8 ■ $— It..”'*1  (77) 

V MiB/  dr  iB  S 4 s dr  ’ 

where  the  coefficients  are  obtained  from  equations  (5b),  (62),  (66),  (70), 
(72),  (7U),  and  (76), 
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HYDRODYNAMIC  DERIVATIVES 


Ihe  theoretical  values  of  the  hydrodynamic  lift  and  drag  force  co- 
efficients and  , as  well  as  their  derivatives,  are  given  in  this 

section.  Tnese  coefficients  are  functions  of  incidence  angle  a , speed 
V , ana  depth  h below  the  smooth  water  surface. 

A hydrofoil  moving  near  a free  surface  will  suffer  a decrease  in 
lift  as  compared  to  one  at  infinite  depth,  due  to  the  reduced  virtual  mass 
associated  with  the  foil  and  the  flow  modification  due  to  wave  formation. 
The  theoretical  expression  for  the  lift  coefficient  of  a hydrofoil  has  been 
determined  by  different  investigators  (e.g.,  References  5 and  6).  Analysis 
and  comparison  of  the  lift  theories  indicate  that  the  one  due  to  Keldysch 
and  Lavrentiev  (Reference  5)  leads  to  results  which  are  mere  consistent 
with  experimental  data. 

The  Keldysch -Lavrentiev  theory  gives  the  following  expression  for  the 
lift  coefficient  cf  a two-dimensional  flat  plate  which  approximates  the 
actual  thin  airfoil: 

CX  « 2nag(Px  - agP2)  , (7b) 


where  a is  the  geometric  setting  relative  to  the  smooth  water  surface 
& 

(two-dimensional  angle  of  attack;,  and  P and  P„  are  dimensionless  ex- 

1 c 2 

pressions  dependent  upon  depth  and  Froude  number,  \/gc  , For  the  small 
angles  usually  experienced  in  the  motions  considered  in  this  report,  the 
two-dimensional  lift  coefficient  may  be  approximated  by 

r..  m p . (?o) 

'gl 

For  an  actual  three-dimensional  foil,  it  has  been  determined  in 
Reference  7 that 


aCL  _ .„r  "A  ♦ 6h/c 

TT  n?l  nA  ♦ (bh/c)  ♦ RnP1 


(bO) 


and  hence, 


■M 


V3S* 
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= 2txcxP-, 


nA  + 6h/c 

nA  + (oh/c)  + ItnP^ 


where  a in  this  case  is  the  actual  angle  of  attack  of  the  foil  relative 
to  the  zero  lift  angle.  This  expression  for  lift  coefficient  can  be  used 
to  determine  the  derivative  of  lift  coefficient  with  respect  to  depth, 


<3Ct/ <3h  , and  thus  can  be  related  to  the  components  making  up  the  longi- 

- — i_J 


tudinal  stability  derivatives. 


The  value  of  dC^/ dh  is  given  as 


L _ 9 I nA  ♦ 8h/c  J 

<5h  n nA  + (8h/c)  + UnP^  dh 


<3P  (32nP1/c)-Un  -j—  (nA+8h/c) 

1 £nA  + (8h/c)  + UnP^j^ 


Since  the  coefficients  are  to  be  evaluated  at  the  equilibrium  position  and 
since  H is  defined  as  the  deviation  from  equilibrium  depth,  the  following 
relation  exists: 


V?F;e  \jrl  ,{WM 


where  the  subscript  "en  denotes  equilibrium  condition. 


Motion  of  the  hydrofoil  below  the  free  surface  results  in  wavemaking 
with  which  there  is  associated  a wave  drag.  This  effect  is  due  to  the 
proximity  of  the  bound  vortices  to  the  surface.  In  a finite-depth  channel, 
there  is  a critical  speed  above  which  there  is  no  wave  drag,  which  is  de- 


fined as 


' VcVffhB 


where  h^  is  the  tank  water  depth.  Since  the  Keldysch-Lavrentiev  theory 
is  derived  on  the  basis  of  infinite  depth,  there  is  no  mathematical  mech- 
anism that  will  make  this  wave  drag  component  disappear;  therefore,  it  is 
necessary  to  determine  whether  the  wave  drag  is  to  be  included  in  an 
analysis  of  stability  tests  in  a finite-depth  channel. 
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The  theoretical  expression  for  the  u&ve  drag  coefficient  is  given 
in  Reference  5 as 


2 2 


” ag 

K 


(«i  - w 


(85) 


where  K is  a form  of  Froude  number,  i.e., 


K 


(66) 


and  and  are  functions  similar  to  the  previous  and  P^  . 

The  same  type  of  approximation  leading  to  equation  (79)  is  used  for  the 
wave  drag  coefficient,  resulting  in 


2 2 

n a 

£.  n 

K 


(97) 


In  addition  to  the  wave  drag,  there  is  the  induced  drag  due  to  the 
trailing  vortices  of  the  finite  span  and  the  profile  drag  due  to  fric- 
tional and  form  drags.  The  expression  for  the  induced  drag  coefficient 
is  given  in  Reference  7 as 


nA  + Bh/c 


(88) 


It  is  assumed  that  the  profile  drag  does  not  vary  much  with  angle  of  at- 
tack or  depth | hence. 


(89) 


The  values  of  profile  drag  to  be  used  in  evaluating  the  total  drag  of  a 
hydrofoil  may  be  determined  from  wind-tunnel  tests  at  the  appropriate 
Reynolds  number. 

The  finul  expression  for  the  drag  coefficient  is  then  the  sum  of 
the  contributions  of  profile  drag,  induced  drag,  and  wave  drag,  i.e., 
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C = C + c + c 
33  DP  Di  °w 


2 2 
n a 


C = C + — t-  + Q 

D Dp  nA  + dh/c  K "1 


With  this  expression  for  , it  is  a simple  matter  to  determine  the  re- 
quired derivatives.  Thus, 


I — — 

(nA  + 8h/c)CL  - hC^/c  n2 a c 

(na  + dh/c) 


K dh 


dC^  hC  -3 — 2n  ci  da 

IT  nA  * dh/c  K Q1  Jt 


Since  a = a - a , where  a is  the  angle  of  zero  lift  which  is  prac- 
g o o 

tically  constant  for  the  range  of  speeds  considered,  <3a  / - 1 and 

g 

hence 

dC 

dC  hC  2n\ 

D = I*  da.  . g n l oV 


The  variations  of  lift  and  drag  coefficients  with  angle  of  attack  and  with 
depth  have  been  theoretically  expressed  and  may  be  evaluated  for  applica- 
tion within  the  expressions  comprising  the  stability  derivatives.  The 
values  of  these  coefficients  and  their  derivatives  for  the  rear  foil  are 
found  in  the  same  manner  as  for  the  forward  foil.  Ary  modification  of 
these  factors  due  to  downwash  is  accounted  for  by  the  presence  of  downwash 
terms  involving  de  /da  . 


ml. 


mt&m 
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DETERMINATION  OF  DOWNWASH  ANGLE  < AND  CF  d</dg 

The  angle  of  attack  at  the  rear  foil  of  a tandem  hydrofoil  system 
is  the  difference  between  the  angular  setting  of  the  rear  foil  (relative 
to  zero  lift  angle)  and  the  mean  downwash  angle  < . Thus, 

*R  " a - e , (9U) 

where  a is  the  difference  between  the  geometric  setting  of  the  rear  foil 
relative  to  the  surface  of  the  undisturbed  fluid  and  the  angle  of  zero  lift. 

The  3en3e  of  t is  positive  when  it  tends  to  reduce  the  angle  of 
attack  at  the  rear  foil.  It  is  detsrmined  from  the  slope  of  the  sinusoidal 
3u*face  wave  generated  by  the  bound  vortex  of  the  forward  foil  (based  on 
two-dimensional  theory).  Beyond  the  first  l/k  wave  length  downstream,  the 
surface  i3  represented  by  a sine  wave  whose  equation  (Reference  10)  is 
given  as 

£ =»  - C ce  e sin  ^ , (95) 

h v‘ 

where  £ is  the  displacement  cf  the  surface,  s is  the  distance  between 
the  quarter-chord  points  of  the  forward  and  rear  foils,  and  h is  the 
depth  of  submersion  of  the  forward  foil.  For  the  small  slopes  considered 
in  the  linearized  wave  theory,  the  downwash  angle  « is  determined  as 
the  negative  of  the  slope  of  the  surface  at  the  particular  point  in  con- 
sideration. The  surface  wave  slope  is  then 


is  found  to 
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SUMMARY  OF  EQUATIONS 


Ihe  differential  equations  of  a disturbed  longitudinal  motion  of 
the  tandem  hydrofoil  system  are  summarized  below.  They  were  derived  as 
equations  (16),  (k2) , (h9),  (61),  and  (77)  and  are  collected  here  for 


convenience. 


3—  • x u ■ x fi  * x 0 - x J 9.  q = 0 

or  u w 0 


- v * v - v - v • 2 # ' l1  ■ 3 


t (V."  »>  £ - s - ^ . .*3  9 . L.^) 4 . 5| $ ■ o 

XB  iE  XB  \ ^iB/  dT  iB 


de  a _ 

~ dr  q 0 


0 - e ♦ - 0 

ar 
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STABILITY  DETERMINANTS 


The  complete  system  of  differential  equations  may  be  solved  by 
expressing  each  of  the  five  variables  as  a sum  of  exponentials  of  the 


power  X_.r  , i.e., 

A A X'T 

A V A X 

U = L U^e 


_ X.r 

P^p.e1 


X.r  >.  t ~\ 

w = X w.  e 1 , 0=Ze.ei  , 


A 

n 

9 '-i 


Z/V  . •• 

A ^ 1 
q . © 


The  differential  equations  are  reduced  to  a system  of  five  algebraic 
equations  for  X , the  values  of  which  determine  the  stability  of  the 
. motion.  In  order  to  have  a consistent  set  of  equations,  the  deter- 
minant of  these  equations  must  vanish.  Now,  by  defining 


Z.  Z1  • 

-i  , "W  w , 
a - 1 + , b 

H-  M 


-ix{ m»  + m' «)  z m'  * 

.b  • j? 


-z 

u 

aX  - z 

w 

-ze 

z ' • 

2^ 

r 

n 

\J 

, •>  ^mw 
bA-^ 

_/^e 

1B 

*b 

m 

dX  -r3- 

0 

0 

-X 

0 

1 

0 

1 

-1 

X 

0 

The  system  is  of  uhe  fifth  order.  The  stability  quintic  will  now 
be  obtained  in  the  usual  way  by  equating  to  zero  the  determinant  of  the 
system  (101).  The  resultant  quintic  expression  may  be  written  as 

X5  + BX1*  + CV*  + Tfk*  EX  + F - 0 , (102 
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where 


am  m z 1 

be  - - z d - w 

iR  w 


ad  + bz  1 •/ it 

q r 


d(z„  + x z - z x ) ♦ olz  - Cx  > 

■ uw  uw  \ 6 u M/ 

ad  + bz 1 »/a 

q 

..  / m-.z'«  m x z'»  \ m 

£ (m  c ♦ am  - -2_a  . .D  q - -2.  (z  + ax  ) 
1B\W  ® /i  ^ / ifi  w u 

2 

ad  ♦ bz ' »/m 

qr 


m.(c  - a)  + z - m z + am  x ♦ m cx  - 
iB  3 w G v 0 0 u w u 

ad  ♦ bz ' •/u t 

q r 


m„x  z'» 
9 u q 


— (msx  -mz-razx  -m  x z ) ♦ b(z  - z x + z x.) 
ig  quw  q 8 wuq  quw  9 0u  u 0 

2 

ad  ♦ bz  1 m/ii 

q 


— f m„z„  - m.z,  * m.z  - m z„  ♦ m z,x  - in  z x ♦ m„x  (a  - c 
1R  L 9 0 69  3 w w 9 w0u  0 w u p u 


ad  ♦ bz'»/fi 

•^■(nux  z ♦ m zQx  - m z xQ)  ♦ b(z  x,  - zQx  ) 
i ' 3 q u q 3 a q u 9 v u 9 P u' 

D 

♦ — 

ad  ♦ bz’ •/«  w 

q 

T£  <Ve  • ViTVV.1  * Tjf  (V»  - Va*  Vi 

ad  ♦ bz  ’ *^/ji  ^ 
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The  method  of  analysis  utilized  to  d-.termine  the  dynamic  stability 
without  actually  solving  the  equations  of  motion  is  by  means  of  the 
Hurvitz  criteria  or  Routh's  discriminant  (Reference  11).  For  the  particu- 
lar case  of  a quintic  expression  in  X , the  following  conditions  must 
hold  in  order  to  have  stable  solutions  (real  part  of  X less  than  zero). 
With  an  expression  of  the  form 


i 3 2 

a^X5  + + + a3^  + ah^  * a5 


the  criteria  of  stability  requires  that 
1)  all  the  ai  > 0 


2>  *1  ao 


*1  ao  0 


*3  *2 


>0  , 


ai  ao 


a3  a2  al  >0  > 


a a a 
3 2 10 


a5  \ a3 


a5  aU  a3  a2 


0 0 


ar  a 

5 b 


In  the  longitudinal  stability  tests  of  Reference  12,  constant  speed 
was  maintained  and  no  longitudinal  surging  motion  was  encountered;  hence, 
the  entire  x-equation  can  be  neglected  for  this  case.  The  remaining  four 
equations  may  then  be  analyzed  in  the  same  manner  as  for  the  three  degrees 
of  freedom  system.  The  stability  determinant  will  be  a quartic  in  this 
case,  derived  frcm  the  determinantal  equation 


aX  - z 


z 

w 

_ze 

-ZB 

3*  X - c 

^®w 

^me 

" s 

s 

m 

dX-jS 

B 

-X 

0 

1 

-l 

K 

0 

(101*) 


Die  quartic  may  be  written  as 
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Xu  * E^a3  ♦ C1K4  ♦ DjA  ♦ 2^  " 0 , 


(105) 


wnere 


B,  = 


am  m z ' • 

be  - —a  - z d - js-j a 

W H-H 


ac  * bz 1 •/ll< 
0 ^ 


z m / mQz ' »\ 

. J * z.b  ♦ — w— a~  - — (m  e ♦ m a - 

® ^3  \ w ® M ^ / 


ad  ♦ bz ' •/ll‘ 

q ^ 


Di " 


(b  ' 1^)  * £ fVe  - V.  " "3(c  - B)] 

_ 

ad  ♦ bz  1 •/u. 

q ^ 


(106) 


E,  = 


(m  z - m z ♦ m.z  - m z.) 
ifi  3 6 9 3 3 w w 0' 


ad  + bz ' ./u  ' 

q r 


The  previously  outlined  stability  analysis  based  upon  the  Hurwilz 
criteria  is  applicable  to  thi3  case  as  well.  When  the  expanded  deter- 
minants are  greater  than  zero,  the  values  of  X will  be  such  that  the 
motion  is  stable  and  the  disturbances  are  damped  out. 
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SOLUTION  OF  EQUATION'S  CT  MOTION 


An  exact  solution  of  the  equations  of  motion  that  will  yield  the 
trajectories  of  the  system  will  be  obtained  through  use  of  the  Laplace 
transform.  The  five  equations  previously  found  can  be  reduced  to  a 
system  of  three  equations  through  elimination  of  the  variables  w and 
q , which  are  related  by  purely  kinematic  conditions  to  the  other  vari- 
ables. Thus,  the  following  equations  are  obtained: 


d-u  _ A i de  . , x de  . 0 


-j — x u - (x  ♦ xje  , 
or  u w 8 dr 


w dr 


a , . , ,d0  2 q d20  - d0  d2B  ~ 

- v < W*  * (a  ’ c)ar  - ^^7'V  a?  - * 0 


(107) 


' i^(mw+ne^  + (b"  i^)af*d 


nq\  d0  . . d28  ^mB  Q ^mw  d 


d?  S K S 


d2S 


dr 


The  solution  for  the  case  of  constant  speed  with  no  longitudinal 
surging  will  be  considered  first.  This  case  leads  to  two  equations  with 
6 and  0 as  variables,  the  entire  x-equation  and  terns  such  as  z^u 
not  being  considered.  The  two  equations  are  then 


-a if-  y*,  af  - * 0 


h (Ws*(b-^)  dJ  - if  ^ « ' b ‘ °; 


(10b) 


Ihe  Laplaoe  transform  will  now  be  applied  to  these  equations.  The  fol- 
lowing relations  define  the  operation  of  the  transform  on  the  variables 
0 and  3 and  on  their  derivatives  (Reference  13): 
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^uiir  i a. ru-> 

TntVr*ma  r>r\ 


r 1 • <T> 

L [e(r)j  = Jo  e " efTjdr  = e(p) 

L ^3 (t  )]  = f* e pr3(r)dr  - 3(d) 

L [|f]  = P0(P)  * 8(0) 


.2 

a e 


dr 


p e(p)  - P0(O)  - 6'(0) 


If]  = p3(p)-  3(0) 


rd2 


rf 

dr 


P 6 ( p ) - P3(0)  - 0* (0) 


where  G(0)  , 0'(O)  , 9(0)  , and  3'(0)  are  the  initial  values  ( ^ “ 0) 
of  0 , de/dT  , 3 , and  d9/dr  , respectively. 


The  transformed  equations  are 


z' • 2 

— £ P - (a  - c.)p  ♦ (z  ♦ z ) 

e(p)  ♦ 

ap2  - y>  + Z3 

M 

i—  J 

- 

3(p) 


z 1 • 

f z • . 

dp(uj  ♦ — 5*  o[0) 

r J 

p ♦ -^6:(0)  - (a  - c)e(O)  - z 3(0)  ♦ a9'(0) 

Lm2  w 

(109) 


dp2  ♦ 


/ m \ 

(b'h)p~h( 


m ♦ m ) 

w e 


e(p)  ♦ 


O Mm 

-bp  ♦ — p - -r-2- 

B B. 


3(p) 


L 


d6(0)  - b3(0) 


(b  - i) 6!0 


fi  m 

) ♦ d6 ' (0,1  3(0)  -bp'(O) 

H 


(no) 


Solving  these  two  simultaneous  linear  equations  for  8(p)  with  the  aid 
of  the  definitions 
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z' . 1 "z> . 

. a?(C)  + — 3 0(0)  p ♦ e:(0)  - (a- c)e(O)  -z  3{0)*a3*(0) 
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2 
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T 
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4? 

p-^) 

-(ad  + bz'*//x2)  (p^  ♦ B-jP3  «•  ^p2  ♦ DjP  ♦ 

where  the  denominator  contains  the  stability  quarti?  previously  deter- 
mined (equation  (105)).  Evaluation  of  the  determinant  in  the  numerator 
then  yields  the  final  value  of  0(p)  i 

e(o)p3  ♦ [b^o)  ♦ e^o)]?^  + fo  + 


p ♦ B-jP"  ♦ Cxp  ♦ D^p  ♦ E1 
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The  functions  0(r)  and  9(f)  can  now  be  evaluated  by  means  of 
the  Mellin  inversion  theorem,  using  an  integration  in  the  complex  plane 
along  the  so-called  Bromwich  path,  with  p ■ u ♦ iv  being  the  complex 

variable i 

u*ioo 


6(t)  “ f e^sfpjdp 

u-i  00 
u*ioo 

9(r)  - ^ / ePr8(p)dp 


u-i  oo 

The  singular  points  of  the  integrand  are  the  zeros  of  the  denominator, 
which  in  this  case  is  the  stability  quartic.  After  determining  the  four 
roots  of  this  quartic  polynomial,  it  is  expressed  as 
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(p  - Xj(p  - X2)(p  - x,)(p  - X,  ) 


wnere  X.  are  the  roots.  The  solution  is  then  given  as 
u*ioo 

e(i 


•^ioo  r , 

1 f dT6^0^-  * |_Bie'°)  + 9'(0)Jp  + 0p 

;t)  'hi  J , * ~TP—  x,5(p  -xj(p  -xj(p  -x. ) 

U-100  J C J U 


<1/ 


-up  (121) 


6(r)  =*  *)  Residues  of  [epTfi(p)j  at  X^ 


i-1 

Similarly,  for  3(f)  , 

h 


3(r)  = > ' Residues  of  [e  PTB(P)]  at  X.^ 


(122) 


(123) 


Thus  the  determination  of  the  trajectories  of  the  motion  ia  reduced 
to  a solution  of  an  algebraic  equation  and  an  evaluation  of  the  stability 
derivatives  based  upon  equilibrium  parameters. 

For  the  case  in  which  there  is  an  additional  degree  of  freedom,  i.e., 
longitudinal  surging,  there  are  three  equations  for  the  three  unknowns  ft  , 
8 , and  6 given  by  equation  (107).  The  transformed  equations  are  then 


fx 

(p  - x^)u(p)  - ft  p * (xw  ♦ x0)  8(p)  ♦ (jyj  - Xg)3(p) 


b 


(0)  - ft  6(0)  ♦ Xft(0) 


) 


(I2h) 


A Z'a  2 

iu(p)  * ftz  P - (a  - c)p  ♦ (\*ze) 


e(p)  ♦ 


ap  - z_^p  ♦ 7. 


3(F) 


z * * 1 

a5(0)  + —ft  5(0)1 p + 


r?  j t 


A * / A \ / . , \ A / a\  n A /7>\  , . A • / A \ 

v • \ v j - - C jvyv ; - 2 o ^ U J ▼ cio  • u y 


(125) 
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r-  / m \ I-  ~ To 

dp^  + ( b - t3.  1 p - f-  (m  + m ) 0(p)  + -bp2  t-r^p  ■ -3  0(p)  = 

L \ J'B  ' XB  W 6 J L aB  XB  J 

“ ~|  H m \ ftm 

d9(G)  - b0(O)  p + (b  - t^]  0(0)  + d9'(0)  ♦ — - 3(0)  - bJ3'(0) 
[\  1B  / XB 


(1^6) 


Since,  initially,  the  surge  velocity  u is  zero,  u(0)  =0  and  equation  (12U)  can  be  re- 


written as 


(p  - xu)£(p)  - (j[f  P ♦ (xw  ♦ xQ)  6(p)  + (x.^p  - Xg)B(p)  = - yt  9(0)  + xw  3(0)]  . (127)  | 


Solving  for  u(p)  gives 


[-  -pt  ®(°)  +xw9(°]  "[J  P + (\  + xe)] 

X |— ^ p2  - (a  - c)p  + ( 


(V  - V 


p - (a- c)p  + (zw  + zfl)  (ap  - zxrfp  + z ) 


dp2  + (b-^P-^(«  +I”J  (-bp^p-^ 

\ S/  iB  ' w °j  \ XB  h) 


(12b) 


(p  - x ) 

U 


p + (x  + X.) 
v ' w e 


<V  - V 


~z'*  2 

—§•  p"  - (a  - c)p  + (zw+  zQ) 


(ap  - z^p  + zj 


r 2 / mo\  U , H / , 2 H-\  W 

I dp  * b - T—  J p - -* — (m  ♦ m,,)  I - op  •>-— p 2. 

L \ V 3B  V eJ  V XB  ’*B  - 


s where  X and  Y aro  given  by  equations  (ill)  and  (112). 


After  expanding  the  denominator,  the  stability  quintic  for  this  case  (equation 

2 

(10 2))  multiplied  by  the  term  -(ad  + bz'^//1-  ) is  obtained.  Evaluation  of  the  mime? 
tor  then  gives 
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bx^z 

dzc  * — ' - - a(x  z - x z ) 
6 fJ.  ' w u u w 


ax  m 

bcx  + — ~ 3- 
u ic 

- 


ad  + bz ' • / uf 

q r 


x m z ' • m z ' • / m \ , 

JLiL-a  + P q _ a (b  - -aU 

MS  w v s/  1 


MS 


. \ 


'B 


m 

w 


6(0) 


ad  + bz 1 */u  ‘ 

q •• 


m z 1 « 

dz  + Sr* 
x + — w-.-.gr, 

u ad  + bz'  ‘/ii1 

q ^ -I 


6'  (0) 


Mam8  hz 

-TT-b2B 


ad  + bz  1 •/a 

q r 


Mamw 

bz  - i 

8(0)  + - £-2  8 ' (0) 

ad  + bz ' »/ii 

q r 


d(xoz 


- x z ) -ft 

u ur  \ 


m \ 

. _a) 


(x  z ■■  x z + z ) 


p u "u  P'  y H U ~U~W  8' 

2 

ad  + bz  1 »/f£. 


•(133) 


x m z 
q w u 


x m z 1 • . . 

u p Q , M(a  ~ c) 


MS 


B 


(rax  + in  ) 
w u p 


6(0) 


ad  + bz 1 •/ll  ‘ 

q 


+ 


x m z 1 * 


u w ^ 

MS 


- d(x  z 

V u 


ad  + 


bz  VM 


X z ) 

u w 


e>  (o) 


(Tu 


) + £ (m0z. 


S 


3 w 


m z + ax  m ) 
w 0 up 


Ri'n't 


ad  + bz'./jH-' 

q 


b(x  z - x z 
w u u w 


V ♦« i <»  * * 

B 


3^’  {°) 


ad  + bz ' */u.  ‘ 

<r 


Security  Information 
CONFIDENTIAL 


confidential 
Security  Information 


u> 


1 


b(x  z - x z ) * 
M u S'  i 


M(a  - c)mp*»qzp 


- (m  x -m  x.)  }e(0) 

h q P g * ) 


ad  + hr. 1 •/ti  ‘ 


d(Vu"  Vp}  '“^r^je,(0) 

'i  ^ 


iB  ( Vg " yj " iB  (Vw " mwz 


ad  ♦ bz 1 •/ li? 

q ^ 


ad  ♦ bz'^/fi 


flax  m ' 

b(x.z  - x z„)  «•  — -u  - 
. ’ » “ 8 *B  j 


8'(0) 


ad  4 bz 1 •/ii> 

q r 


“'W^W 

2 

ad  + bz 1 •/ii 

q r 


6(0)  d(a  - c)  4 ?(‘-4)  8'(0) 


ad  4 bz 1 •/u ' 

q ^ 


z m bz  x 

rf_ 

lB 


bz  4 -.y  3 4 - dx  z — ♦ m - mQ) 

6 iD  A1  w u 1 i ' w 0 P' 


ad  ♦ b*'«/u 

q r 


B X 

2 


ba  P(0) 


be 


am 

a 

S 


, ad 


bz » •/a 

q n 


3'(o) 


s>  b<°) 
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f[xu(\  * V - 2U<\  * V]  * -f/  (aw  * ”e> 


ad  ♦ bz'  •/u‘ 

qr 


^(a  - °>t\  ♦»,)  -(■»  - (=tt  * »„>} 


ad  + bz  1 • /u^ 

q r 


x z' . / 


ad  ♦ bz ' ./a 

q r 


x m z x m z 


ad  ♦ bz 1 */m 

q 


ax  (m  ♦ m.)  - (a  - c)x  m„  ♦ z (m  ♦ in  ) 

uv  w 6'  ' U 8 wv  w e' 


ad  ♦ bz  '•/u 


A ♦ b(z  . ♦ z J 


ad  ♦ bz ' •/a 
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fX  2 / jjj  \ 

-4--  (m  ♦ m ) + [z  (x  ♦ x ) - x (z  ■*  z )~j  ( b — ] 

iE  v e'  Lu  »»  e u'  w e'J  \ ±BJ 


ad  ♦ bz ' •/ 

q ^ 


1 

— — (a  - e)(mw  ♦ nQ)  V 0(0) 


ad  ♦ bz1 ;/u 

q ^ 


/d[x  (z  ♦ z ) - 2 {a  x )~j  + _H__a  (m  + m ) \ © * I 

|_u  w 6 u v 6 J fi i ' v B j 


(132) 

cont  • 


ad  + bz  1 -/u‘ 


tr~  fz  (x  m - m xn)  ♦ X m(z  * zA)  - x z (m  + ni  )"1  3(0) 
iB  lU  W 6 W 6 u 3 w 6 u w w 6 J 


ad  * bz ' . /u  ‘ 

q 


fbfz  (x  > x ) - x (z  ♦ zjl  ■*•  (m  ♦ m ) \ 3'  (0) 

L Lu  w e u w e'J  i ' w e j 1 ' 


ad  ♦ bz?  ,/il 

q r 


The  solutions  lor  u(r)  , e(r)  , and  3(t)  are  found  in  exactly  the  same 
manner  as  previously  for  the  case  of  two  degrees  of  freedom.  These  solutions 
are  given  by 


«w  ■ E Residues  of  ['  Tu(p)]  at  \± 
i=l 

0(r)  ™ y Residues  of  [yjr9(p)J  at  X 


(131*) 


(135) 


? _ _| 

0(r)  -E  Residues  of  |^e^r3(p)j  at  X^ 


(136) 


where  X.,  are  the  solutions  of  the  stability  quintic  for  this  case  (equation  (102)). 


mm  * 
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ATTLICATICN  OF  EQUATIONS  OF  MOTION  TO  LONGITUDINAL  STABILITY  TESTS 


The  previously  derived  equations  of  motion  (99)  will  now  be  applied 
to  actual  experimental  studies.  The  longitudinal  stability  tests  of 
Reference  12  will  be  used.  Since,  in  these  tests,  constant  speed  was 
maintained,  the  whole  x-equation  can  be  loft  out  of  the  calculations. 

The  resulting  equations  for  this  case  are  given  by  equation  (10b)  and  the 
stability  determinant  is  given  by  equations  (10l),  (105),  and  (106). 

The  tests  were  made  on  a tandem  system  made  up  of  two  aspect  ratio  20 
hydrofoils  of  rectangular  plan form  with  0°  dihedral.  The  physical  dimen- 
sions and  conditions  that  remained  constant  in  the  tests  are  listed  below: 


X=  § = U.167  ft. 


W = 63.3  lb. 


2.5  

C = - U.CUOJ  IT,. 


13.51 

qF,R  " 12 


1.12b  ft. 


S = 1.736  ft.? 

S_  - SD  - 0.86ti  ft.2 
F R 


w 

fJL=  — a — = O.iliOU 
gp-SZ 


The  equilibrium  depth  was  one  chord,  hg  = 0.2083  ft.,  and  equilibrium  trim 
was  0°,  Equilibrium  speed  was  in  the  neighborhood  of  8.50  ft. /sec.,  but 
was  slightly  different  in  each  of  the  cases  tested.  The  moment  of  inertia 
of  the  cys  tciiij  the  forward  and  rear  pre-set  angles  of  attacx,  and  the  loca- 
tion of  the  longitudinal  center  of  gravity  were  also  varied  for  each  case. 
The  center- of -gravity  location  was  determined  by  the  equilibrium  conditions 
necessary  for  a static  balance,  equations  (l)  and  (2). 

The  first  ease  to  be  treated  is  that  of  the  far-forward  center-of- 

gravity  position,  x./£,  = 0.25  . Equilibrium  speed  V was  8.10  ft. /sec., 

-*•  2 2 e 

the  moment  of  inertia  was  such  that  i„  = K_  / JL  = 0.3601  . and  t.hp  angles 

d a • ' ° 

of  attack  (relative  to  zero  lift)  were  a„  - 0.1715  radian  and 

F 

ae^  = 0 radian  . He  hydrudynamic  coefficients  and  derivatives  are  to  be 
determined  on  the  basis  of  the  above  conditions  and  will  then  be  used  to 
determine  the  stability  derivatives. 
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For  this  case,  the  stability  quartic  is  determined  and  the  coeffi- 
cients of  the  powers  of  X are  evaluated  with  the  following  results: 

a0  = 1 a^  = Di  = O-llifl; 

a,  - B = 5.5093  a,  = E - 0.0505  . 

1 hi 

a2  ' C1  “ 631*5 


Application  of  the  Hurwitz  criteria  shows  that  all  the  a^>0: 


“ a^a^  - aQa  = 25.1865  > 0 


*i 


ao  0 


al(a2a3  " alV  " a0a3£  = 7 * 1919  > 0 


Since  the  determinants  are  greater  than  zero,  the  values  of  > are  such 
that  the  motion  i3  stable  and  the  disturbances  are  dajrped  out. 

The  roots  of  the  quartic  are  found  to  be 


X,  - -0.91*8  X3  - -0.0335  ♦ 0.101*3  I 

X2  - -h.h9h  X^  - -0.0335  - 0.101*3  i . 


The  existence  of  the  conjugate  complex  roots  indicates  oscillatory  sta- 
bility, which  was  discovered  to  be  the  condition  in  the  stability  experi- 
ments for  this  configuration. 

For  this  case,  the  initial  values  of  the  variables  were  a positive 
trim  of  2°  and  a submergence  of  two  chords,  so  that  the  hydrofoil  system 
was  initially  displaced  from  equilibrium.  The  initial  values  of  the 
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variables  and  their  derivatives  were: 


0(0)  = 2°  = 0.031*9  radian  ; 


= 0'(O)  - 0 


8(0)  = -0.3561  ; 


. 


dr 


T =0 

= fl'(O)  = 0.031*9 


r=0 


By  solving  the  equations  of  motion,  the  trajectories  are  round  to  he 
6(r)  = 0.0037  e“,5>il8r  - 0.0001*2  e~^"^T 

+ e“*033^T  (0.03393  cos  0.101*3^  + 0.0091*  sin  0.101*3  r)  (13?) 


P(r)  = 0.05261  e~‘9lxtiT  - 0.01615  e"^,li9hT 

- e“°033^r(o .1*5991*  cos  0.101*3 r - 0.39662  sin  0.10l*3r)  . (138) 

These  solutions  are  plotted  in  Figures  1 and  2 and  indicate  a low  ampli- 
tude, low  frequency,  danpea  oscillation,  which  is  in  agreement  with  the 
experimental  findings.  They  are  close  to  the  trajectories  which  would  re- 
sult from  an  exact  solution  of  the  equations  but  differ  slightly  because 
of  the  small  numbers  dealt  with  and  the  manipulation  of  the  complex  numbers 
to  determine  the  coefficients  of  the  exponential  terms.  This  can  be  seen 
in  Figures  1 and  2 at  r = 0 where  the  initial  conditions  are  not  exactly 
satisfied. 

The  next  case  to  be  investigated  here  is  the  one  with  the  center  of 

gravity  located  at  Xj ^/JL  0.569  . The  value  of  V was  8.1*5  ft. /sec., 

the  value  of  i^  was  0.211*8,  and  the  angles  of  incidence  were 

aQ  “Qp  “ 0.101*7  for  the  forward  and  rear  foils.  Evaluation  of  the 
eF  eR 

stability  determinant  for  this  condition  gives  the  following  values  for 


= B1  = 5.1*51*5 

ig  ■ = 6*- 9506 


a3 

% 


D1  “ 0.728? 


E1 


0.0201 
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Tnt’^  values  then  satisfy  the  first  requisite  for  stability  u.-Jei  the 

Hurwit2  criteria,  i.e.,  that  all  i.  > 0: 

1 


a. 

1 


a1a2  - aQa  - 37.1633  > 0 


ao  ° 


0 ah  a 


= a^a^  - a^)  - aQa  « 26.U977  > 0. 


Since  the  determinants  are  greater  than  zero,  the  notion  in  this  case  is 
stable. 

by  assuming  the  same  initial  conditions  for  this  configuration  as 
for  the  previous  case  of  the  far-forward  center  of  gravity,  and  solving 
the  equations  of  motion,  the  trajectories  are  found  to  be 

0(t)  - 0.01602  e~*°^T  «.  0.02005  e',06tiT 


- 0.00115  e'1* 765 T-  0.00001  e~3'^7T 


(139) 


S(t)  = 1.9707  e"'°^r-  2.35535  e~*06dT’ 

♦ 0.03W1  e"1,7&5T  - 0.00658  e"3,557T  (D*0) 

plots  of  which  are  given  in  Figures  3 and  u. 

For  the  case  of  the  center  of  gravity  located  at  x^/i.  ■ 0.60U  , 

V was  6.37  ft. /sec.,  i^,  was  0.216b,  and  the  angles  of  incidence  were 
e d 

ae  - 0.0673  radian  and  ae  - 0.1222  radian.  The  stability  quartlc  is 

F ^ 

doterreined  and  the  coefficients  are 
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a^  = Dx  ■=  0.6107 

a = B1  « 5.5020 

a,  = E = 0.0029 
4 1 

a2  - C1  - 6,0801* 

These  values  of  a^  > 0 satisfy  the  first  requisite  for  stability  under 
the  Hurwitz  criteria.  The  required  determinants  will  also  be  greater  than 
zero,  indicating  stability  for  this  case.  However,  the  value  of  the  coef- 
ficient E^  is  of  such  small  magnitude  that  in  the  region  \mder  considera- 
tion, there  cannot  be  any  definite  statement  regarding  the  stability  of 
the  motion,  because  of  the  definition  of  as  the  sum  cf  products  of 

small  numbers.  The  actual  motion  of  this  configuration  in  the  stability 
tests  was  a slow,  unstable,  diving  divergence  and  was  very  close  to  the 
divergence  boundary,  the  actual  location  of  which  was  approximately  at 
*jJJL  ~ 0.60. 

In  order  to  shew  the  trend  of  the  solutions  of  the  equations  in 
agreeing  with  the  experimental  results,  a static  balance  was  found  by- 
using  the  data  of  Reference  ll*  at  a center-of-gravd ty  location  given  by 
V*  - 0 .735  * The  equilibrium  speed  was  6.50  ft, /sec.  and  the  moment 
of  inertia  was  chosen  such  that  i ~ 0.1(091*  . The  angles  of  incidence 

D 

on  the  forward  and  rear  foils  were  a„  = G.O696  radian  and 

eF 

ae  * 0?1571  radian.  The  stability  quartic  is  evaluated  for  this  case  and 
the  coefficients  are  found  to  be 

aQ  = 1 a^  = D_l  = 0.1*169 

al  = B1  = ^#38  al*  “ E1  = "°'°133  • 

a2  - ci  - 3.7112 

Since  a^  is  a negative  number,  the  first  requirement  of  the  criteria  for 
stability  is  not  satisfied  and  hence  the  motion  is  unstable. 

The  effect  of  an  additional  degree  of  freedom  (surge  motion)  is  now 
investigated  for  the  stable  case  at  = 0.569  * The  resulting  sta- 

bility quintic  has  the  following  coefficients: 
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a = jl)  = 1.0667 

j 

a,  * K = 0.0336 
a 

o s K n '1  O vn 

w5  * • 

Valid  conclusions  as  to  the  stability  of  this  configuration  are  not  possible 
because  of  the  small  order  of  magnitude  of  the  coefficient  F . However, 
the  fact  that  F is  very  small  is  somewhat  indicati/e  of  a lessening  of  sta 
bility  due  to  this  additional  degree  01  freedom. 


= B = 5.';7b3 
a,,  = = 7.17*9 
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CONCLUSIONS 

The  equations  derived  in  this  report,  which  are  based  upon  the  as- 
sumptions outlined  on  page  6 , adequately  describe  the  motion  of  a tandem 
hydrofoil  system  in  smooth  water.  Analysis  of  the  resultant  equations, 
with  the  coefficients  evaluated  on  the  basis  of  theoretical  hydrodynamic 
formulas,  indicates  a.  duplication  of  the  type  of  motion  experienced  in 
tank  tests.  The  trend  in  the  motion  with  oscillatory  stability  for  a 
far-forward  center-of -gravity  location,  darrped  stable  motion  at  a center- 
of-gravity  position  near  the  middle  of  the  configuration,  and  divergent 
instability  at  points  further  back  indicates  the  validity  of  the  theory 
developed  herein.  Therefore,  for  constant  speed,  the  derived  equations 
are  sufficient  for  investigating  the  longitudinal  dynamic  stability. 

A preliminary  study  on  the  effect  of  including  the  x-equation  in 
the  analysis  for  the  case  of  varying  speed  indicates  a lessening  of  sta- 
bility. This  additional  degree  of  freedom  in  longitudinal  surging  is, 
however,  not  too  important  in  the  present  case,  i.e.,  motion  under  an 
undisturbed  free  surface.  It  is  of  great  importance  though  for  the  case 
of  notion  in  waves  where  the  waves  cause  variations  ~'n  the  forward  speed. 

The  equations  developed  in  this  report  can  be  extended  to  a study 
of  motion  in  waves  and  controlled  motion,  by  employing  the  proper  types 
of  forcing  functions  on  the  right-hand  side.  Solution  of  the  resulting 
equations  in  this  case  by  use  of  the  Laplace  transform  method  will  in- 
volve only  additional  algebraic  operations.  Further  investigation  of  the 
present  system  for  varied  equilibrium  conditions  will  enable  an  investi- 
gator to  determine  the  proper  stable  range  of  operation  for  a particular 
design  having  fully  submerged  hydrofoils  in  tandem. 
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PITCH  ANGLE  CHANGE  VS.  TIME  (NONDIMENSIONAL) 
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